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Abstract 

The renormalization group equations for the general 2x2, complex, neutrino 
mass matrix are shown to have exact, analytic solutions. Simple formulas are given 
for the physical mixing angle, complex phase and mass ratio in terms of their initial 
values and the energy scales. We also establish a (complex) renormalization invari- 
ant relating these parameters. The qualitative features of the physical parameters' 
renormalization are clearly illustrated in vector field plots. In both the SM and 
MSSM, maximal mixing is a saddle point. 
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1 Introduction 



Recent experiments [|l], Q have revealed some salient features of the neutrino mass matrix. 
It is now rather well-established that the neutrino masses are tiny, and that at least 
some of the mixing angles are large, or even maximal. Considerable efforts have been 
devoted to a theoretical understanding of these features. Any theoretical model faces an 
important issue. Is the model valid at the high or low energy scale? How are the physical 
parameters related between the two scales? These questions are answered by studying the 
renormalization group equation (RGE). In addition, RGE may offer a natural mechanism 
which can drive the neutrino mixing angle to its maximal value. 

The RGEs governing the neutrino mass matrix were worked out some time ago 
and have been studied extensively (see e.g. @, ||). Most studies considered the case of 
real matrices, although some dealt with complex ones also. In this paper, we wish to 
point out that there is an exact, analytic solution for the RGE of the general 2x2 mass 
matrix. Simple formulas are given for the physical mixing angle, complex phase and mass 
ratio in terms of their initial values at an energy scale. We further show that there is a 
(complex) RGE invariant relating the running of these three variables. The fixed points 
of the physical parameters are obtained and their stability is determined. The nature of 
the running near these fixed points is clearly illustrated in vector field plots [0. 

In both the standard model (SM) and the minimal supersymmetric standard model 
(MSSM), the effective Majorana neutrino mass matrix can arise from a dimension-five 
operator. Its RGE (in the basis where the charged-lepton mass matrix is diagonal) is 
given by (t = In fi/M x ), 

-^ m v = —(Km v + m u P + P T m v ) , (1) 

where k depends on a combination of the coupling constants and, for two flavors (for 
definiteness, we consider and v T ) ||, after absorbing the muon Yukawa coupling term 
into K, 

P = P T = x (l- a 3 ) , (2) 

f (l£-I©/4, SM 
X -\ _(#-jg)/2, MSSM W 

Here y T = \/2m T /v and y T = y/2m T /(v cos/?) are the r Yukawa couplings in the SM and 
MSSM respectively, with v ~ 246 GEV and tan/3 is given by the ratio of the two Higgs 
VEVs in MSSM. Also, y^ and y^ are similarly defined. 



2 Solution 

A formal solution || to Eq. (P is 



m v (t) = e - K V <T3 m l ,(0)e ? ' J3 (4) 
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where we have ignored the i-dependence of the coupling constants^ so that / ndt ~ Kt, 
etc., and 

k' = k + 2 X (5) 
£ = Xt- (6) 

It is convenient to factor out the determinant 

m u = v /m 1 m 2 M , (7) 

then detM = +1, and the mixing angle and the complex mass ratio are contained solely 
in M. We have 

= e-"'Vmi(0)m 2 (0) (8) 
M(t) = e ?ff3 M(0)e ?<T3 . (9) 

The overall scale ^Jm^rn^ has a simple dependence on t. The mixing angle and the 
(complex) mass ratio evolve via a transformation that is similar to the effects of the 
familiar seesaw model 0. The difference lies only in that while |£| ^> 1 for the seesaw 
model, for the RGE usually |£| < 1 (~ for the SM). Nevertheless, the general 

analytic solution obtained earlier [§, |] is equally valid here. We will now discuss this 
solution in the context of RGE. 

As was shown earlier ||, a symmetric and complex 2x2 matrix with det = 1 can be 
parametrized in a standard way. We write 

M(t) = U(t)e~ 2vaa U(t) T (10) 

U(t) = e -i^ e -i9^2 e -i^ 3 ^ ^-Q 

where r\ is given in terms of the mass eigenvalues by 

,= ih,=*. (12) 

4 m 1 

At t — 0, M(0) and U (0) are similarly defined in terms of (770, #0, «o, 0o)- Here mi and m\ 
are, by definition, positive definite. The relative phase of the mass eigenvalues is given 
by 40 so that M becomes real for = and <p = n/4, corresponding to the same sign 
and opposite sign mass values, respectively; 9 is the physical mixing angle. The phase a 
can be absorbed into the arbitrary phase of the charged leptons and is not observable. 
However, as emphasized before (see refs. @] and |§), it evolves with t and so influences 
the evolution of the other parameters. We can rewrite M = M(t) as |J 

_/ e- 2ia (ch.2fj - c 2e sh.2fj) -s 2 e sh2r/ \ . . 

I -s 2e sh2f) e 2ia (ch2fj + c 2e sh2^) ) 1 ' 



1 To include the i-dependence, one needs only to replace Kt by J Kdt etc, in the appropriate formulae 
in the following. Numerically, we have verified that the difference is negligible. 
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where we use the notation ch2r7 = cosh2?7, c 2 # = cos 26, etc. and r\ — rj + icj). 

The RGE evolution of the parameters is given by M(0) — > M(t). From Eq. (P), it is 
obvious that the off-diagonal elements of M do not evolve. We have thus the (complex) 
RGE invariant: 

s 2e sh2fj = s 26o sh2?7o (14) 

where the subscript denotes the values at the high energy scale t = 0. We note that in 
terms of the physical neutrino masses, sh2f/ = (^Jm^/rnie 12 ^ — ^/mi/m 2 e~ l2(?i ) /2. 

Eq. (§) also implies two complex relations between the diagonal elements of M(t) and 
M(0). They can be used to solve for the four unknowns (a, 6, <ft, 77) in terms of their initial 
values (a , 6 , (f> , r] ). After a short calculation we find 

tan2£ = W(c 2 A Q ch2Q (15) 

c 2 e — S fi tanh24 + S/tan2Aa 

tan2A« = ^ (16) 

_ v , . v 1 - (mi/m 2 )g - 2i(m 1 /m 2 ) o g 40o 
C ° th2 "° = Ss + lS '= l + ( mi /m a )8-2 (mi /m a W (17) 

where Aa = a — ocq, the S's are real, and the masses in the last equation are at t = 0. 
These solutions agree with Ref. |J, where they were obtained by a somewhat different 
method. We do not write down the solution for fj, since, knowing 6, 77 and <fi can be 
determined by Eq. (|T4"D. 



3 RGE 

The RGEs for the physical parameters can be worked out. From Eq. @, the RGE for 
M(t) is 

j t M = X {M,a 3 }. (18) 
We may decompose M into Pauli matrices, 

M = T^Mm , (19) 
where the coefficients Mj are complex, and M 2 = since M is symmetric. Then 

± Ml = (20) 

^(M„ ± M s ) = ±2 X (M„ ± M3) . (21) 

Eq. (p0|), which says that the off-diagonal elements of M do not run, is just another form 
ofEq. (0), 

j t (s 2e sh2f]) = . (22) 
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The diagonal elements of M run according to 

^[e- 2iQ (ch2r/-c 2e sh2fj)] = +2 X e~ 2ia (ch2fj - c 2e sh2??) (23) 

^[e 2ia (ch2fj + c 2e sh2fj)] = -2 X e 2ia (ch2f] + c 2e sh2fj) . (24) 



Eqs. fl23|) and ( p4j ) can be readily rewritten in terms of RGE for the various parameters. 
It is found that, for the unobservable phase a, 

d = J^ (25) 
dt shAr] v ; 



The RGEs for the physical parameters are 

dt] 
~dt 



\ c 26 (26) 



d< P S H 

= - XC29 s^ (27) 

dt 



X -529 [c 2 coth 2r/ + s 2(j) tanh 2r/] . (28) 

It is easy to verify that Eq. fl22|) is satisfied. Also, the last equation agrees with the known 
result for the case of real mass matrix |3], f|, |5| 

d9 m 2 + m x 

-77 = XS29 (29) 

at m 2 — mi 

valid for same sign (<fi = 0) and opposite sign (<f) = tt /4, or mi — > —mi) mass val- 
ues. Our result, with the convention of positive definite mi and m 2 , interpolates these 
two equations. Note also that the equations can be expressed in terms of tanh27y, with 
[dL coth 2r\ + s 2( p tanh 2rf\ = 1/Re(tanh 2fj) and s^/shArj = Im(tanh2f/)/Re(tanh2f/). Fi- 
nally, the right-hand sides of Eqs. (|25| - |2"8"D are independent of the unobservable phase a, 
as it should be. Explicitly, from Eq. (f|), e^°" 3 commutes with e iacr3 , so that any RGE 
evolution does not depend on a. 

It is interesting to note that, for = 0, there is a simple geometrical interpretation 
for Eqs. ( p6|) and fl28|). As was shown earlier, by writing M(t + dt) = e^ s e~ 2riS ' n2e e^ a3 
with £ = xdt, Eq. (|9[) can be regarded as a velocity addition problem (8v + v + Sv) with 
v = (tanh2r/)n 2 6i and 6v = —xdtn 3 . Here h 2 e = c 2 en 3 + s 2 ehi. The resultant Lorentz 
boost has the rapidity 2(r/ + dr]), with dr\ = 5v\\ = ~xdtc 2 g, while its angle with respect 
to the third axis is increased by 2d6 = 2Sv±/\v\ = 2xdts 2 g/ tanh 2rj. 



4 Phase Portraits 

To better understand the qualitative effects of the RGE, we rewrite Eqs. (p6|-|28|) in terms 
of a slightly different combination of the physical mass eigenvalues || 

+ uo- m 2 e 2i ^ + mie' 2 ^ 

Z = COtll27? = — r-r (30) 

m e 2«t> _ m e -2i<t> 
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In terms of this variable, the RGEs are 

dz . {-l+z 2 )z* 

It = 2x °»7h (32) 

where z* is the complex conjugate of z. 

The neutrino mass matrix, Eq. (|i~3"D, possesses various symmetries. For example, under 
the transformation 

- 1 - e (33) 

accompanied by an interchange of the mass eigenvalues 

f] — > —f] (m2 <-> mi) (34) 

0^-0 (35) 

or, equivalently, 

z -> (36) 

the diagonal elements are invariant while the off-diagonal elements change sign. However 
the sign of the off-diagonal mass-matrix element, which is the RGE invariant (Eq. ([14]) ), 
is not a physical observable since it may be absorbed in the unphysical phases with a 
redefinition of the neutrino wave function (u^, v T ) —>■ (u^,—u T ). Thus the mass matrix 
and the RGE evolution equations respect this symmetry. 

The evolution equations in terms of complex z and 9, Eqs. (|3l|) and (|32|) , have five 
fixed points. The evolution is always towards or away from one of these fixed points, 
i.e. the parameters do not evolve towards infinity. Of the five fixed points, only three are 
physically distinct in that the above mass matrix symmetry maps two fixed points into 
two others. We have obtained the stability of the fixed points for t decreasing (i.e. running 
towards the infrared) in the SM by finding the eigenvalues of the Jacobian (see e.g. ||). 
The attractive fixed points are 9 = 0, z = +1 and 9 = 7r/2, z — — 1 which corresponds to 
no mixing and a massless muon- neutrino. The repulsive fixed points are 9 = 0,z = — 1 
and 9 = tt/2,z = +1 which corresponds to no mixing and a massless tau-neutrino. The 
final fixed point is a saddle point, attractive in some directions and repulsive in others, and 
it is at 9 = 7r/4, z = 0. This point corresponds to maximal mixing with equal magnitude 
but opposite sign mass eigenvalues. These stabilities are for the SM, for the MSSM all 
evolution directions are reversed, so the attractors and repulsors are reversed also. 

These RGEs and the fixed points are graphically displayed in Figs. |I]-|3| where the 
direction fields are plotted for t decreasing. Starting from some initial point specified by 
the high energy theory, the evolution of z and 9 follow a trajectory in these spaces and 
the plotted unit arrows show the directions tangent to this trajectory. The direction field 
is independent of t and \. The trajectories can be found by taking the RGE equations, 
Eqs. ( |3~1~D and (|32|), dividing one by the other and integrating to get 

^ - C (37) 
z z — 1 
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where £ is a constant. This equation is just the square of the RGE invariant, Eq. ([b 

If the mass matrix is initially real, then RGE evolution preserves this. The evolution 
of real parameters, lm(z) = 0, is shown in Fig. [I]. The symmetries of the neutrino mass 
matrix, Eqs. ( |33"D andfl36|), are apparent. The dark curves in Fig. [I] correspond to £ = — 1 
in Eq. fl3"7p. These curves connect all of the fixed points. Note that half of these dark 
curves flow towards the maximal mixing fixed point and end there, the only trajectory 
that does so. Trajectories with ( < — 1 or £ > are those on the left and right of the 
plot and they can pass through maximal mixing. Trajectories with — 1 < ( < connect 
positive and negative z values and never attain maximal mixing. 

The general behavior for lm(z)^ is shown in Fig. [|, where the vector field is plotted 
in the full, 3-D parameter space. This figure is for the range < Re(2;), < Im(2;), 
< 8 < 7i / 2, so the 'floor' of this figure corresponds to the right half of Fig. |l|. The 
behavior in other regions of the parameter space can be obtained from this by symmetry 
operations. The three fixed points are shown using a dark sphere, light sphere and a cuboid 
to denote the attractor, repellor and saddle point. Note that the figure resembles the 
vector field of a dipole, with evolution away from the repellor and towards the attractor. 

The evolution equations are singular at Ke(z) = if lm(z)^ 0. Thus large changes in 
the mixing angle are possible then, despite the small size of x- The behavior in this region 
is shown in Fig. |3|. Here the direction field is plotted for the lm(z), Re(z) plane. Note 
that this direction field is independent of the magnitude of 8, except that the direction 
of flow does depend on the sign of COS26 1 . The fixed points for 8 = and 8 = 7r/4 are 
included in the plot. As Fig. |3] shows, the flow across the Re(z) = value does not occur if 
Im(z) 7^ 0. The flow on the right-hand side is towards the stable fixed point there. On the 
left-hand side, the flow is away from the repulsor, but at the same time the mixing angle is 
increasing and eventually exceeds 8 = tt/A at which point the flow direction reverses, the 
repulsor becomes an attractor and the evolution approaches the now stable fixed point 
by retracing its path in the lm(z), Ke(z) plane. 



5 Analysis 

We turn now to a quantitative analysis of our results. First, we consider the case of real 
mass matrices, 0o — or 7r/4, so that Aa = from Eq. ( |16| ) and 

tan 28= S ^° h \ - (38) 
C29 ~ tanh 2£ 

where, from Eq. ( |TTD 

(mt+mx) = cot h277 , O = 



tanh2?7o , 4>o = tt/4 . 
£| < e>(l(T 2 ) for the A 

ately that for O = 7r/4 (opposite sign mass eigenvalues), 8 ~ 8 . For 4> Q = (same sign 



Since |f | < 1, (|f | ~ 1(T 5 for the SM and |£| < O(l0~ 2 ) for the MSSM), we see immedi 
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masses), tan 2$ exhibits the resonance behavior, and maximal mixing is obtained if 

c 29o = tanh2£/tanh2?7o , (0o = 0) . (40) 

For a generic value of 9 < ir/4, this condition can be fulfilled provided that £ and rjo are 
of the same sign, and that rjo is very small (nearly degenerate masses). In terms of the 
masses, if we define 

6 = 1 - (m 1 /m 2 ) (41) 
then the resonance condition becomes 

c Wo ~ 4£/5 (0 O = 0) . (42) 

For SM, £ < 0, it is seen that 9 < 6 , if m 2 > mi. That is, under "normal" situations, 
RGE drives the mixing angle toward zero as the energy decreases. For MSSM, £ > 0, we 
see that 9 can become large if there is near degeneracy, m 2 — mi (m 2 > mi). Of course, 
if 6*o > 7r/4 (which amounts to the interchange mi <-> m 2 ), or if mi > m 2 , the direction of 
the running of the mixing angle is reversed. These conclusions agree with earlier results 
obtained by numerical integrations of the RGE|. 

At the same time, the RGE invariant (Eq. ([14])) gives the running of the mass ratio 

s 2e sh2r? = s 2 o sh2??o , (0o = 0) (43) 
s 2e ch2ri = s 2do ch2^ , (0o = 7r/4) . (44) 

While rj ~ r] for O = 7r/4, for O = and at the resonance, s 2 q ^> s 28o , it is seen that 
sh2r/ <C sh2r/o, *-e., the masses are driven toward even more degeneracy as 9 becomes 
maximal from a small 9q. 

For complex mass matrices (</>o 7^ or 7r/4), we need to use Eq. (|i~5|). From Eq. (f2§|), 
it is seen that for x ^ lj d9/dt is appreciable only if 77 <§C 1. In other words, RGE effect 
can be large only if 5 = 1 — [m\jm 2 ) 1, as for real matrices. We will confine our 
discussion in this parameter region 5 <C 1. Further, for definiteness, we assume that £ > 
and 7/0 > 0. (The other possibilities can be easily analyzed.) For finite values of 0o, it is 
easy to see that 

E R ~ S /(2s 2 2 J (45) 
£/ ~ -cot 20o (46) 

where higher order terms in 5q are neglected. The condition for maximal mixing is then 

^ - (c 2 e 4<j> + 4JC26 ) , (00 ± 0, tt/4) (47) 

This reduces to Eq. ( ^2| ) as O — > 0. However, when O = 7r/4, we had found earlier 
that there is no resonance. This means that the limit 0o — > 7r/4 is actually very delicate. 

2 Note also that, for general values of cf>, similar behaviors can be deduced from Eq. ( p8| ) 
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As we will see, what happens is that the width of the resonance also shrinks to zero as 
0o ~~ > 7r /4- This behavior is borne out in the numerical calculation shown in Figs. ^ and 

Although the analytic formulae for 9, r\ and </> have been explicitly given, as we have 
seen in the previous analyses, their detailed behaviors are quite intriguing. It is useful 
to have an overview of these functions. For this purpose we present, in Figs. |] and [5], 
numerical calculations of the function 9(t). Traditionally the plot of 9(t) is made with 
t as the variable. However, as shown in Eqs. ( (4~2"D and fl4"?|), the interesting region is 
determined by £o/£ ~ ^(1)- So an equivalent way is to plot 9 vs. 5o (or (mi/m 2 )o), while 
keeping t as a constant parameter. In Fig. [|, we present a 3D plot of sin 2 29 vs. (mi/m2)o 
and 40o- We take a positive £ = \t = 1CT 2 (corresponding to MSSM with tan/3 = 40 
and t = j^j\n(m z /M x ) } M x ~3x 10 10 GeV). This sets the scale in {mi/m^o, with 
(1 — mi/rri2)o = So ~ 10~ 2 . Thus for a given initial value 9o = n/12, the plot gives sin 2 29 
as a function of the initial parameters (mi/ 1712)0 and its phase 40o- The resonant behavior 
of sin 29 is obvious. The position of the maximum is, for 4> = 0, at 5 ~ 8£/V3. The 
position of the maximum shifts with changing <f) , and is well described by Eq. ((47|). 

Fig. |5] shows 2D slices of the 3D plot in Fig. |. We plot, for various values of 4> , 
sin 2 29 vs. (mi/m 2 )o- It shows clearly the shifting positions of the maxima as 0o changes. 
Also notice the change of width as a function of 0o- It can be shown that, provided 
4£ 2 -C c 2 ^ tan 2 20 o , the width is given by 

A(m 1 /m 2 ) ~ 8£c 2 ^ (W4j\/ c 24>o + > ( 48 ) 

so that the resonance disappears as 4>o — * tt/4, in agreement with Eqs. ( p8|) and (|3^). 
Finally, the singular nature of the limit <fio —> tt/4 is confirmed when we compare the 
curves with O = 0.97r/4 and O = 7r/4. 

We do not show 3D plots for rj and 0. Since from the RGE invariant, 

s 2 8 sh2?7 c 2 = s 2 8 sh2r] c 2( j> , (49) 
s 2 o ch2?7 s 2( p = s 2 e ch2?7o s 2<j}0 , (50) 

once we know 9, they can be easily determined. Note that for nearly degenerate masses 
(i] ^ 1), it can be seen from Eq. ( [50D that there is an anti-correlation between 9 and 
with the phase <fi being minimal at maximal mixing. Also, at maximal mixing and with 
nearly degenerate masses, simple expressions for the neutrino mass ratio and the relative 
phase can be deduced from Eqs. ( f4"9"D and (pO]), 

S24, — s 2 e Q S2<t> (51) 
5 ~ cW,^/^ 1 - 4e 4<f> > (52) 

where 6 = 1 — mi/mi- As can be seen from Eq. (|52]), for small #0, the neutrino masses 
become more degenerate at maximal mixing. 
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6 Conclusions 



We have examined the RGEs for the physical, two-flavor, neutrino parameters: mixing 
angle, mass ratio and their relative phase in the SM and MSSM. These equations turn out 
to have relatively simple forms which we analyze in detail. The qualitative nature of the 
evolution is clearly illustrated in our vector field plots. A more quantitative description 
is also given using our exact, analytical solution for the evolution. In addition, we have 
also found a complex RGE invariant which correlates the running of the three physical 
parameters. 

The phase portraits show the direction of the evolution throughout the parameter 
space. The fixed points and their stability are included on the plots. It should be noted 
that maximal mixing (6 = 7r/4) is a saddle point in both the SM and MSSM, and the 
infrared stable point (attractor) in the SM (MSSM) corresponds to no mixing and a 
massless muon-neutrino (tau- neutrino). Thus to get maximal mixing at the experimental 
scale requires a particular choice of mass ratio and phase at the high-energy scale. 

The parameter choices that give maximal mixing have been calculated with our exact 
solution. As is well-known for real parameters, large evolution requires nearly degenerate 
masses. With the addition of a complex phase, the peak position shifts and the resonance 
region becomes narrower, so to achieve maximal mixing at the experimental scale requires 
a finer tuning of parameters at the high-energy scale. 

In the literature, considerable interest has been focused on the possibility of generating 
large mixing through RGE running. The detailed solution for RGE suggests both oppor- 
tunities and limitations for this scenario. Since the end result depends very sensitively on 
the initial conditions, any such model must be treated carefully. In particular, one needs 
to study the three flavor problem in detail. Although it is very difficult to obtain exact 
solutions in this case, approximate solutions seem well within reach. We hope to return 
to this problem in the future. 
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Decreasing Energy, Standard Model 




Figure 1: Phase portrait for the Standard Model RG equations when lm(z) = 0. The 
different fixed points are shown with solid circles, open circles and grey square denoting 
attractors, repellors and saddle point. The solid curve shows the trajectories that connect 
the fixed points. The arrows are reversed for the MSSM. 

ff 




Figure 2: The SM RGE evolution in Re(z), 9, Im(z) space, for the range < Re (2), 
< Im(z), < 9 < 71/ 2. The different fixed points are all confined to the lm(z) = 
plane. They are shown using a dark sphere, light sphere and a cuboid which denote the 
attractor, repellor and saddle point. 
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Decreasing Energy, Standard Model, Cos (29) > 
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Figure 3: The SM RGE evolution in the Re(z), Im(z) plane for cos 2$ > 0. The = 
fixed points are shown with solid circle, open circle and grey square denoting attractor, 
repellor and saddle point. For cos 29 < 0, the figure is similar except all directions (and 
stabilities) are reversed. Note that 9 also evolves. 
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Figure 4: 3D plot for sin 2 29 vs. initial values (mi/m2)o and 40 o . The inputs are fixed 
at 6 = tt/12 and the RGE factor f = 0.01 (MSSM). 
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Figure 5: sin 2 29 vs. (m 1 /m 2 ) with inputs 9 = vr/12 and f = 0.01 (MSSM). The solid, 
dotted, dashed, and dot-dashed curves correspond to 40 o = 0, n/2, 0.97T and it respectively. 
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